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We show that the local in-gap Greens function of a band insulator Go(e, ky, rj_ = 0), with ri the 
position perpendicular to a codimension-1 or -2 impurity, reveals the topological nature of the phase. 
For a topological insulator, the eigenvalues of this Greens function attain zeros in the gap, whereas 
for a trivial insulator the eigenvalues remain nonzero. This topological classification is related to 
the existence of in-gap bound states along codimension-1 and -2 impurities. Whereas codimension-1 
impurities can be viewed as ’soft edges’, the result for codimension-2 impurities is nontrivial and 
allows for a direct experimental measurement of the topological nature of 2d insulators. 
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The topological characterization of condensed states 
of matter has emerged as a prominent research inter¬ 
est over the last few decades. The flourishing of the 
Quantum Hall effect (QHE) [1] in particular elucidated 
many connections between physical signatures and topo¬ 
logical invariants [ 2 ], which supplement the order param¬ 
eters of the usual symmetry-breaking Landau-Ginzburg 
paradigm. More recently, however, it became apparent 
that topological order can also arise by virtue of symme¬ 
try, and in particular the very common and robust Time 
Reversal (TR) symmetry is sufficient to establish the ex¬ 
istence and stability of topological insulators [3-5]. This 
is quantified via a Z 2 invariant, and results in gapless he¬ 
lical edge states or chiral Dirac fermions localized at the 
perimeter of the sample, in two and three dimensions, 
respectively. The topological insulator has proven ex¬ 
tremely rich both experimentally and theoretically [9, 10]. 
The concept has been generalized to a “periodic table” 
describing various discrete symmetries and dimensions 
[ 6 - 8 ]. Lattice symmetries can similarly lead to further 
topological distinctions, result in “crystalline topological 
insulators” [ 11 ], for which a general classification has been 
provided [ 12 ]. 

One may ask whether the topology of band insula¬ 
tors has some local signature? In fact, in this Letter, 
we will show that even the fully local in-gap Greens 
function contains information about the band topology, 
which is then directly accessible by experiments. The 
natural way this insight arises is through the study of 
impurities[13-15], similar to how the space group classi¬ 
fication can be probed using lattice defects[16-23]. Con¬ 
sider a codimension -1 impurity line or surface in an insu¬ 
lator. In the limit where the impurity strength diverges, 
V —> 00 , such an impurity acts like a real edge which, 
following the bulk-boundary correspondence, should host 
zero gap metallic bound states in the topological phase. 
For finite V the codimension-1 impurity surface can thus 
be viewed as a ” soft edge”. The codinrension-2 impurity 
lines or points do not host gapless states in the strong V 
limit, so a priori there is no reason to expect they probe 


topology. However, we will see that they in fact inherit 
the topological structure of the ’’soft edges”. This non¬ 
trivial result implies that by probing bound states around 
a point impurity in a two-dimensional insulator, one can 
experimentally distinguish between the topological and 
trivial phase. 

Mathematically, the theory for impurity bound states 
shows that the existence of bound states in the topolog¬ 
ical phase is directly related to zero eigenvalues of the 
local in-gap Greens function. Consequently, we propose 
that the presence or absence of zero eigenvalues in the 
local in-gap Greens function is a signature of the band 
topology. 

The remainder of this Letter is organized as follows. 
We first introduce the model and the theory of impurity 
bound states. Then we show that the local in-gap Greens 
function, relevant for codimension -1 and -2 impurities, 
has zero eigenvalues if and only if the phase is topologi¬ 
cal. We then propose an experiment that directly probes 
the Greens function zeroes for a point impurity in a 2d 
insulator. Finally, we relate our results to the known Z 2 
classification and the bulk-boundary correspondence. 

Theory of impurity bound states- We begin with a 
translationally invariant system, described by a mini¬ 
mal time-reversal invariant two-band model. The generic 
Hamiltonian assumes the form 

H o = J2 c la d n(. k )Ta0 c ^0 (!) 

k otft 

where are the 4x4 Dirac gamma matrices satisfying a 
Clifford algebra. We choose 7 0 = cto< 8 >T 3 and 7 * = <Ti®T\. 
Here, the a and r Pauli matrices act in the spin and or¬ 
bital space, respectively. Time reversal (TR) symmetry 
then implies that do(k) must be even and dj(k) must 
be an odd function. In particular, we focus on the rep¬ 
resentative cases that do(k) = M — 2B JT(1 — coski ) 
and <ij(k) = sinfc.; [24-26]. This is the familiar class of 
models displaying topological nontrivial regimes for pa¬ 
rameter range 0 < M/B < 4d. This specific choice is not 
expected to restrict our results, as a topological insulator 
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generically has by adiabatic continuity the form of a lat¬ 
tice regularized massive Dirac Hamiltonian. Moreover, 
extra terms that respect TR symmetry will not change 
the results described below. 

The real frequency Greens function in the gap reads 


G 0 (w,k) 


1 

w - d M (k)y^ 


L0I4 + d A 1 (k) 7 M 
w 2 ^|d(k)| 2 ' 1 ’ 


Subsequently, we introduce an impurity into the system, 
which in general can be described by the Hamiltonian 


H V = Y1 c L V <xfl( r ) c rfl- (3) 

ra/3 


To find the corresponding spectrum in the presence of the 
impurity, one needs to solve the (differential) Schrodinger 
equation 


(H 0 + H v )ip t ( r) = eip t ( r) (4) 

where e is the energy of the state. This can be trans¬ 
formed into an integral equation, [27] 

^e( r ) = ' 52 G 0 (e,r-r')V{r')ip e (r'). ( 5 ) 

r' 

For an insulator, the real-frequency Greens function is 
well-behaved inside the gap and displays exponential de¬ 
cay as a function of r. This implies that an in-gap so¬ 
lution of Eqn. (5) will yield a bound state around the 
impurity. 

The existence of an in-gap bound state depends on Hq 
and the shape of the impurity potential V(r). However, 
the qualitative difference between topological and trivial 
band insulators is found in Greens function features, and 
is therefore largely independent of the choice of impurity 
potential. Let us therefore consider the simplest possible 
choice: a constant V (r) along a n-dimensional plane in 
a d-dimensional system (hence codimension d — n). The 
d-dimensional position vector r can be split into the per¬ 
pendicular coordinates rj_ and the parallel coordinates 
r || , so that the impurity potential is given by 

V(r) = V o d ",.* 0 ( 6 ) 


However, an impurity that couples to the density, gener¬ 
ally does not give rise to a diagonal expression, 

Hy = J d d rV(r)^(r)tp(r) (7) 

= ^2 c L c i0 J d d rV(r)0?Jr)&0(r) 

ia/3 

The shape of V 0 can be restricted, though, using symme¬ 
try principles. For example, when we consider nonmag¬ 
netic impurities, TR invariance applies to the impurity 
potential as well [33]. As a result, the matrix Vo has 
only six degrees of freedom, 

V 0 = VI + V 0 a° ® t 3 + Via* 0 r 2 + V 4 a° <g> r 11 . (8) 

Recall that in this notation, time reversal is T = ia 2 K 
where K is complex conjugation. Parity, on the other 
hand, is given by 7 0 = a 0 (g> r 3 . If we require both parity 
and TR the form of Vo is even further constrained to 

V 0 = Vl + Vo7°. (9) 


Since the translational symmetry is not broken along 
directions parallel to the impurity, the impurity bound 
states have a well-defined parallel momenta k||. Thus 
oc e lk ii' r n and the integral equation Eqn. (5) re¬ 
duces to an eigenvalue equation for each ky, 

det[G 0 (e,k||,r x = 0)V o - l] = 0. (10) 


We immediately notice that for the case V 0 = VI, the 
existence of bound states is directly related to the eigen¬ 
values of the local (r_L = 0) in-gap Greens function 
G 0 (e, ky, rj_ = 0). 

Codimension-1 impurities- Let us now consider 
codimension-1 impurities, that is a surface in d = 3 and 
a line in d = 2, having only one perpendicular direc¬ 
tion rj_ = 0, see Eqn. (10). The corresponding Greens 
function in the gap, integrated over the perpendicular 
momentum, can be decomposed in terms of 


9n( e ’ k ||) 
5(e,k||) 


f dk _l d M (k||, k±) 

J 2tt e 2 — |d(k||,fc_i_)| 2 ’ 

f dk± e 

J 2it t 2 - |d(k||, A_l)| 2 ’ 


( 11 ) 

( 12 ) 


where we have used a Kronecker delta to signify our use 
of lattice models and introduced the 4x4 Hermitian ma¬ 
trix V 0 . We note that even if the potential V(r) just 
couples directly to the electron density, V 0 is not nec¬ 
essarily diagonal when expressed in terms of the second 
quantized operators. Explicitly, ip(v) is Ike electron held 
and can be expanded in orbital wave functions as ip( r ) = 
XOi Q c ia0ia(r), where 4>i a { r) is the wave function of the 
a-orbital at lattice site i satisfying the orthonormality 
condition J d d r(/)* a (r)(j)jp(r ) = 5ij5 a p. Hence, the total 
density is diagonal, N = f d d r^(r)ip{r) = J2ia c l a c ia- 


so that G 0 (e, k||,rj_ = 0) = g(e, ky)l + g M (e, ky^. 

At any TR-symmetric point for the parallel momen¬ 
tum, for example ky = 0 or 7 r, the g\\ are vanishing. 
Additionally, g± vanishes since the integrand is an odd 
function of k^. At TR-symmetric points we thus only 
need to consider g 0 and g. We note that this still holds 
if we add e.g. Rashba spin orbit coupling terms, that 
are odd functions of the momentum, to the bare Hamil¬ 
tonian. This allows to us to verify the following results 
also in the absence of any other symmetry but TR sym¬ 
metry. 
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M/B= -0.3 M/B= 0 M/B= 2 

lm x 0 Im x Q Im x 0 



Ml B= 4 M/B= 4.1 M/B= 6 

Im xq Im xq Im xq 



FIG. 1: The flow of the poles of Eq. (15) as function of 
M/B. The two negative residue poles are indicated with 
a dot, whereas the stars mark the poles having a positive 
residue. In the trivial system, two poles with either positive 
or negative residue are located inside the red unit circle. For 
M/B — 0, 4 there are two poles with opposite residue located 
on the unit circle, signaling the transition. The topological 
phase is characterized by two poles of opposite residue in the 
unit circle. 


We will now show that the eigenvalues of the local in¬ 
gap Greens function G 0 (e, kj|, r± = 0), 

A±(e) =g(e)±g 0 (e), (13) 


have the shape displayed in Fig. 2. The difference be¬ 
tween a topological insulator and a trivial insulator, is 
whether the in-gap Greens function has a zero eigenvalue 
or not. Consequently, this result implies that for each 
choice of impurity potential Vo, a topological insulator 
will have at least one in-gap bound state. 

Note that the denominator in the integrand, e 2 — |d| 2 , 
is always negative. In the trivial phase do(k) does not 
change sign throughout the Brillouin zone. This implies 
that go(e) does not change sign. Additionally, because 
e < Mo(k)|, we have g(e) + |^o(e)| > 0 and g(e) - |g 0 (e)| < 
0. Therefore, in the trivial phase, the Greens function 
G 0 (e, k||,rj_ = 0) never has an eigenvalue equal to zero 
for all momenta ky and energies e. 

On the other hand, in the topological phase, when k| 
is chosen as a TR-symmetric point kg associated with 
the projection describing the topological phase [ 12 ], the 
Greens function satisfies G 0 (e = 0,k| = kg,r_L = 0) = 
0. To prove this, we evaluate 


So( 0 , kg) 


r dk dp(fc,kg) 

J-n 2 tt \d 2 (k, kg)| 

r dk M(kg) — 25(1 — cosfc) 

J-n sin 2 k + (M(kg) — 213(1 — cos k )) 2 ’ 


where M( kg) = M — 2 B (1 — cos(fcj)) in terms of 

the coordinates kg. Substituting x = e lk , the integral 
becomes a contour integral over the unit circle and we 


obtain 


1 / dx B + (M ~ 2B)x + Bx 2 

2th t =1 X f ( x ) 


(14) 


with f(x) = (| + 6 B 2 - ABM + M 2 ) x 2 + 2 B(M - 

2B)x(x 2 + 1) + (B 2 — j) (x 4 + 1). This expression is 
then solved by application of the Cauchy residue theo¬ 
rem. The poles are located at, see Fig. 1, 


x 0 


2 B-M± \/l - ABM + M 2 
2B±1 


(15) 


When the system becomes gapless for M = 0 or M = AB , 
we find that there are two zero’s located on the unit circle. 
Thus, from the analytic structure of the integral we infer 
that the regime 0 < M < AB is topologically distinct 
from M < 0 and AB < M. Focusing on the topological 
phase 0 < M < AB , we observe that the residues of the 

, 2B-M±Jl-ABM+M 2 , , , . . , ,, 

two zeroes Xo,± = - 2 B+ 1 -located inside the 

unit circle cancel, 


Res(a; 0 ,+) + Res(a? 0 ,-) = 0 . (16) 

Therefore, in the topological phase, go(e = 0) = 0. When 
M < 0 or AB < M, that is in the trivial phase, the inte¬ 
gral never equates to zero. Together with the universal 
divergence of the both g(e) and go(e) at the band edges 
as shown in the Appendix, we arrive at the generic de¬ 
scription as shown in Fig. 2. In the topological phase 
the Greens function switches sign and thus has at least 
one energy e for which it has a zero eigenvalue. 

Consequently, for any impurity strength a topolog¬ 
ical insulator will always have in-gap states along a 
codimension -1 impurity, whereas for a trivial insulator it 
depends on specific details of the impurity and the insula¬ 
tor. The codimension-1 impurity can thus be understood 
as a ’soft edge’. 

Codimension-2 impurities- The above results on the 
structure of the eigenvalues of G 0 (e, ky, rj_ = 0 ) in the 
codimension -1 case can be extended to codimension -2 im¬ 
purities. In this case, there are two perpendicular direc¬ 
tions kj^ = (Ad[,£r[), 


M e >k||) 

S(e,k y ) 


dk^dk^f 


(27T) 2 

dk^dk^ 


(27T) 2 


d»( k|j,k.i_) 

- M(k||, kj_)| 2 ’ 


e 

- M(k||, kj _)| 2 ’ 


(17) 

(18) 


for p = 0,1,2, 3. It is clear that for any of the perpen¬ 
dicular directions Q± = 0 , as the integrand is odd. 

For the trivial phase, we can show that the eigenvalues 
of Go(e, k||,rj_ = 0) are nonzero throughout the gap, 
since the two-dimensional integral can be done by first 
integrating in one direction, which yields the results from 
the codimension -2 impurities, and then integrating along 
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FIG. 2: The eigenvalues of the Greens function 

Go(e,fc|| =0,rx=0) in the two and three dimensional 
M/B model, relevant for codimension-1 defects. In the 
trivial system (dashed lines, M/B = —1) the eigenvalues A± 
are nonzero, which means for strong enough impurities there 
are no in-gap bound state possible. In the topological system 
(solid line, M/B = 1) the eigenvalues cross zero energy at 
t = 0 and hence in-gap bound states always exist. 


the second direction. Therefore, G 0 (e, k||,rj_ = 0) is 
never zero in the gap. 

Of more interest is the question of existence of zero en¬ 
ergy eigenvalues in the topological regime. Let us focus 
on the 2 -dimensional case, so that there are no parallel di¬ 
rections: we are directly probing the local, on-site Greens 
function. We expect that the terms Q(e) and Go(e) will 
diverge close to the band-edge. In fact, these divergences 
are captured by expanding around the point where the 
gap is minimal, k t ;, 

| d(k x , k v ) | 2 = A 2 + a{k x - k x G ) 2 + b(k v - k v G f + ... (19) 


The diverging part of the integral is then captured by the 
integral 


f dk x dk y 1 

J (27t) 2 2 A Se + a(k x ) 2 + b(k y ) 2 


f° + dq q 

Jo 2'K\[ab 2A5e + q 2 


log Se 
47t\/ ab 


( 20 ) 

( 21 ) 


Hence, G(e) ~ and So(e) ~ in prox¬ 

imity of the valence band. The dependence of the gap 
A on do(k) proves that both eigenvalues in the topo¬ 
logical phase diverge to positive infinity at the valence 
band edge, and to minus infinity at the conduction band 
edge. Consequently, in the topological phase the Greens 
function eigenvalues must be zero somewhere in the gap. 
Details are provided in the Appendix. This proves that 
in d = 2, the completely local in-gap Greens function 
G 0 (e,r = 0) has zero eigenvalues if and only if the sys¬ 
tem is in the topological phase, see Fig. 3. 

This result carries over to the case of line impurities 
in d = 3 topological insulators, if the remaining parallel 
momentum is chosen at one of the TR-symmetric points. 


FIG. 3: The in-gap eigenvalues of the local Greens func¬ 
tion Go(e, r = 0) in the two dimensional M/B model. In 
the trivial system (dashed lines, M/B = — 1) the eigenvalues 
are nonzero, whereas for the topological system (solid line, 
M/B = 1) each eigenvalue is zero for some energy. 



FIG. 4: Typical energy of impurity bound states in a two- 
dimensional insulator as a function of impurity strength V. 
Here, we discern between the topological regime (solid lines, 
M/B = 1) and the trivial regime (dashed line, M/B = — 1). 
For strong V, bound state in the trivial phase has disappeared 
in the conduction band, whereas the bound states in the topo¬ 
logical phase remain. 


Experiment- The existence of these zero eigenvalues 
can be probed directly in experiments, using impurity 
bound states as solutions to Eqn. (10). Imagine a two- 
dimensional insulator, where at one isolated point a tun¬ 
able gate voltage is applied, serving as the impurity po¬ 
tential V. Then using tunneling spectroscopy, the pos¬ 
sible bound states around this impurity can be found. 
Upon increasing the impurity potential V, the energies 
of the bound states shift: for a trivial insulator, one can 
make a bound state disappear into one of the bands by a 
sufficiently strong potential. However, our results show 
that for a topological insulator, for all strong V there will 
always be two bound states. Explicitly, the energy of the 
bound state as a function of V is shown in Fig. 4. 

Classification and the bulk-boundary correspondence- 
The odd number of crossings per spin branch of 
Go(e, k||, rj_ = 0) with the zero eigenvalue axis is a topo- 
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logical property [28-30]. It reflects that in the nontrivial 
regime the system has an odd number of Kramers degen¬ 
erate edge states on either side of the surface and hence 
may be regarded as a consequence of the bulk-boundary 
correspondence [9, 10]. In particular, the bulk TR Z 2 
invariant is in this case simply the product II r . & of the 
parity 70 eigenalues C over the TR points in the Bril- 
louin zone [31]. Moreover, the two relevant poles have 
the same residues but multiplied by sign of the mass, i.e. 
the parity 70 eigenvalues. Hence, only if the choice ks 
in the projected plane is associated with two masses of 
opposite sign, meaning that this cut features a band in¬ 
version, the poles cancel in the above integral rendering 
a zero eigenvalue. This is in accordance with the space 
group classification [ 12 ]. 

For example, for d = 2 the model Eqn. (1) exhibits 
a T (T — pAmm) phase for 0 < M/B < 4 and a M 
(T — pA) phase for 4 < M/B < 8 . From the above con¬ 
siderations we find that M has to satisfy 0 < M < 4 
for G 0 (e,k||,r^ = 0) to develop zero eigenvalues. Tak¬ 
ing subsequently projections onto k x and k y and using 
that M( k s ) = M — 2B JY(1 — cos(fcj)), we thus conclude 
that in the T phase these kg choices correspond to an 
inversion at k = (0,0), whereas in the M phase the in¬ 
version is at k = ( 7 r, 7 r). We stress that this analysis still 
holds if we add Rashba terms, that are odd functions 
of momentum. Similarly, we may add a next nearest 
neighbor term d 0 (k) = —B[ 1 — cos(fc x ) cos(fc y )], d x (k) = 
cos (k x ) sin (k y ), d y ( k) = — sin (k x ) cos (k y ) to the Hamilto¬ 
nian [12], allowing for an additional X — Y (p4) topolog¬ 
ical crystalline phase [11]. This phase is associated with 
the inversion momenta k = ( 7 T, 0 ) and k = ( 0 , 7 r). An 
identical calculation then shows that indeed projections 
alongs k x and k y yield zero eigenvalues of G 0 (e, ky, rj_ = 
0) for both ks = 0 and n. These ideas carry over di¬ 
rectly to three dimensions. Consider for example a T 
(T — pm3m ) phase with an inversion at k = (0, 0,0) for 
0 < M/B < 4, in the projected plane one should now 
choose kg = (0,0) or on a line ks = 0. 

Conclusions and Outlook- We have shown that topo¬ 
logical band insulators can be characterized by the ex¬ 
istence of zero eigenvalues in the local in-gap Greens 
function Go(e,k||,rj_ = 0), where rj_ is the position 
vector perpendicular to a codinrension -1 or -2 impurity. 
Whereas the codimension-1 impurities can be viewed as 
soft edges, the nontrivial result for codimension -2 impuri¬ 
ties suggests one can experimentally probe the difference 
between a topological insulator and a trivial insulator 
using a tunable localized impurity. 

We made some simplifying assumptions in the proof 
presented above, but the results are robust. For example, 
adding more bands to the system, further away from the 
Fermi level, might introduce extra impurity bound states 
whose energies depend strongly on the impurity strength, 
but one can show that these do not generally remove the 
states arising from the low energy bands. Furthermore, 


we showed (e.g. Fig. 4) the persistence of bound states 
for strong potential V but neglected Vo; however, while 
the shape of bound state energy versus potential strength 
changes if V 0 is included, the conclusion that they per¬ 
sist is independent of the ratio Vq/V. Furthermore, the 
principle of adiabatic continuation suggests our proposed 
classification applies equally well to gapped interacting 
systems with a quasiparticle description[32] or topologi¬ 
cal superconductors [ 10 ]. 
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Appendix 

Details codimension-1 case 

We have shown in the main text that Go(e = 0, k|| = 
ks,ry = 0) = 0 only acquires zero eigenvalues in the 
topological regime. Here, we elaborate on the generic 
structure which results in the schematic representation 
as conveyed in Fig. 2. Here we address the behavior 
of G 0 close to the band edges. First, note that g(e) 
is a decreasing odd function in e, and go(e) is an even 
function in e. Close to the band edges, both g(e) and 
go(e) will always diverge as a square root, g ~ 1/y/Se. 
Whenever the gap is located at a T-symmetric point, 
where d,(k) = 0, the strengths of the divergences in g(e) 
and go(e) are equal. Without loss of generality, to prove 
this we choose rn < \ where the gap is located at k = 0. 
We expand |d(k|| = 0, k±)\ 2 around that point, 


|d(k| t = 0, fc_L)| 2 = to 2 + (1 - 2m) fc 2 + .... (22) 

The divergence of g(e) near the valence band edge can be 
isolated by integrating only over a small region (—a, a) 
around the top of the valence band for small Se = e+ |m|, 


g(e = - \m\ + Se) 


f a dk ± \m\ 

J_ a 2tt 2|m|Je + (l-2m)kf' ’ 


m 


8(1 — 2 m)Se 


O(Se 0 


(24) 


A similar argument applies to g 0 ( e ) close to the valence 
band, as the magnitude of the gap |do(k|| = 0, k± = 
0)| = |m|. Therefore ga(e) and g(e) have exactly the 
same divergent behavior close to the band edge. Hence, 
g + \go | only diverges close to the valence band and g — \go\ 
diverges close to the conduction band. Note that this 
analysis is valid for rn < , when the band gap is at 

k = (0,0), and for m > -y, when the band gap is at k = 
("7T, 7r). This implies that this cancellation of divergences 
is present in the trivial phase (which supports the notion 
that there the eigenvalues do not change sign), but also 
for some region in the topological phase. When m > 1/2 
and m < 15/2, the gap is not located at a T-symmetric 
point but rather at (0 ,kc) or (7r, kc). In this case, the 
gap satisfies A 2 = do(^G) 2 +sin 2 kc > do(ko) 2 . This last 
one is important, because then at any T-symmetric point 
of k||, the eigenvalues of G 0 (e, k|) will both diverge to 
positive infinity at the valence band edge. To see that, 
we expand 


|cZ(k)| 2 = A 2 + a(k± - kc) 2 + ■ ■■ 


(25) 
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SO that the divergent parts close to the valence band edge Magnetic codimension-1 impurities 

of the Greens function terms are 


9(e = 

— A + <5e) ~ 

|A| 

(26) 

i/8a| A | <5e 
-d 0 (0, fe G ) 

3o(e = 

— A + Se) ~ 

(27) 

i/8a| A|<5e 


Because |A| > Idol, we find that divergences do not cancel 
and both eigenvalues of G 0 (e, ky) will diverge to positive 
infinity at the valence band edge. This means that at the 
conduction band edge, both diverge to negative infinity, 
therefore implying an energy at which G 0 (e) has zero 
eigenvalues. 

Finally, note that the eigenvalues of G 0 (e, k|,rj_ = 
0 )Vq, relevant for the impurity problem, are given by 


A gv = V 0 g 0 +Vg± 


\ 


(Vg 0 + V 0 g) 2 + )(g 2 -g 2 o) 


(28) 

and are thus directly dependent on the eigenvalues of 
G o (e,kf,rx=0). 


Let us shorty address the fate of codimension-1 mag¬ 
netic impurities. In this case, TR invariance does not 
constrain the specific structure of the H-matrix. Let us, 
however, focus on the simplest case: V 0 ~ (j i 0T°. We as¬ 
sume the impurity line orientation to be the 2-direction, 
so that in d = 2 we have 


G 0 (e, k|j,r_L = 0) = g(u)l + ffoMy 0 + 32 M 7 2 . (29) 


In addition, we assume the following impurity shape 

2 

Vo = ^I4(T i 0T 0 . (30) 

i=1 

which differentiates between the parallel Vi and perpen¬ 
dicular Vi magnetization directions. The resulting struc¬ 
ture has some nice analytic properties. In particular, it 
follows from 


A av(e) = ± y/(g 2 + gl)(V? + V?) + g 2 2 {V? - If) ± 2^/jg ^V 2 + V 2 ) - g\ V?)(gZ(V? + V 2 2 ) + gpf), 


that for every eigenvalue A Gv(e) there is also an eigen¬ 
value —Acy(e)- Therefore, if Acy(e) = 1 for a given e, 
then there is also a solution at — e. We thus conclude that 
impurity bands resulting from magnetic impurities are 
symmetric in energy. This is corroborated with explicit 
tight-binding results in final section of this Appendix. 

In the special case where the magnetic impurity is 
aligned to the impurity, effectively setting Vi to zero, 
the eigenvalues simplify as 

Acv(e) = (±g ± \Jgl+g 2 )V 2 (31) 

which is just ±1 times the solution for the identity matrix 
impurity, Vo ~ 1. Hence, in these cases one just retrieves 
a doubling of the results found in the case of nonmagnetic 
impurities. 

Details codimension-2 impurities 

We here corroborate the ideas of the main text with 
the specific arguments, showing that G 0 (e, ky, rj, = 0) 
in two dimensions has zero eigenvalues in the topological 
regime. From Eq. (20) , it follows that Q(e) ~ 

and Qo(e) ~ ^ close to the valence band. If the 
gap is not at a T-symmetric point, this automatically 
implies that the divergences do not cancel and we are 


left with eigenvalues that all diverge at both band edges, 
leading to the fact that the Greens functions eigenvalues 
have to be zero somewhere in the gap. 

On the other hand, if the gap is at a symmetric point 
we need an extra argument to show the result of Fig. 3, 
hence consider 0 < m < Trivially, Q(e) > 0 for e < 0 
and Q(e = 0) = 0. At the same time, we know that 
both Go{t) an d G(e) have a logarithmic divergence at the 
valence band. However, the difference is that because 
do(k = 0) > 0, the function Go(e) goes to minus infinity. 
At the same time, G(e) goes to plus infinity because in 
that case the numerator is e < 0. Thus the two lines 
must cross if G o(e = 0) > 0, and consequently there is a 
point where G 0 has zero eigenvalues. 

It remains to prove that, in this construction, Go if- = 
0) > 0 for 0 < to < 2 • Tins actually straightforward. 
Since Go{e = 0) = f -^y3o(e = 0, fey) we directly infer 
that 


/ elk 1 

ar KV^I 2 >0 

(32) 

for all fey given to < 2, and hence Go( e = 0) > 0 in the 
desired region 0 < to < 



Tight-binding results 

The main consequence of the difference in structure of 
G 0 (e,k||,rj_ = 0) in the trivial and topological regime 
is reflected in the formation of bound states in the gap 
between the two relevant topologically active bands. As 
this provides for a derived experimental signature, let 
us therefore corroborate the results with explicit tight- 
binding calculations. Fig. 5 shows the results for the two 
dimensional M — B model in presence of a codimension- 
1 impurity given by the identity matrix with strength 
—6 in units of B. The formed impurity bands in the 
gap of the valence and conduction band are clearly vis¬ 
ible. When the impurity strengths are subsequently in¬ 
creased, to V = —100, only the topological regime still 
hosts impurity bands as demonstrated in Fig. 6. This is 


in accordance with the results in the main text. 

The inclusion of Raslrba type couplings to the Hamil¬ 
tonian that do not close the band gap, 

-^Rashba ~ ^ ^ Rp y ( 1 T z ) ((7y sin k x (7 x Sm/Cy), (33) 

k 

results in the lifting of the degeneracy of the impurity 
band as shown. Nonetheless, the qualitative results are 
the same. Some representative cases connecting with the 
previous results are summarized in Fig. 7. 

Finally, the last Fig. 8 shows the spectrum in the 
presence of magnetic impurities. As indicated by the 
above analysis one finds in this case symmetric impurity 
bands. Also in this case Rashba terms do not affect the 
outlined analysis. 
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FIG. 5: In gap impurity bands for various values of M/B for a impurity line that hosts a delta potential having strength 
Vo /B = —6 x 1. Specifically, these simulations were performed on a two dimensional 60 x 90 system with periodic boundary 
conditions. 
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FIG. 6: Analogous plot as the previous one, but now for Vo /B = —100 x 1. We observe that in the gap there are only impurity 
bands for the topological regime. In addition, a difference between the T phase (0 < M/B < 4) and M phase (4 < M/B < 8) 
can be discerned. Namely, in the first case the bands are related to the inversion of the bulk bands at k = 0 and the bands are 
symmetric around the project k = 0 value. Similarly, the M phase is associated with k = r). 
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M/B= 4.8 M/B= 4.8 



FIG. 7: Spectra in the presence of a Rashba term, that break the particle hole symmetry of the original Hamiltonian. The 
strength of the Rashba coupling equates in all cases to Ro/B = 0.6. The resultant effect is the lifting of the degeneracies of the 
impurity bands. Nonetheless, the qualitative features relate to the outlined results and previous figures. Specifically, the panels 
on the left show the system in presence of an impurity potential Vq /B = —6 x 1 , whereas the right have Vq/B = —100 x 1 . 
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M/B= 4.8 M/B= 4.8 



FIG. 8: Impurity bands originating from the inclusion of a codimension-1 magnetic impurity. In all cases the impurity strength 
V/B equated to —6. The left panels show the magnetic result in absence of Rashba terms.The right panels display the same 
information in presence of a Rashba terms with coupling strength Rq/B = 0.6. We observe that the bands are symmetric 
around zero energy. Interestingly, the degeneracy of the nodes is lifted in the trivial regime but not in topological regime. The 
instances shown are for the case of a parallel magnetic impurity, i.e Vo = Via 1 ® r° with i the direction of the line impurity. 
The results for perpendicular realizations resulted in the exact spectra. 



































